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We present a deep optical spectrum of TNJ0924-2201, the most distant radio galaxy at z = 5.19, obtained with FOC AS on the Subaru 
Telescope. We successfully detect, for the first time, the C iV/U549 emission line from the narrow-line region (NLR). In addition to the 
emission-line fluxes of Lya and C iv, we set upper limits on the N v and He n emissions. We use these line detections and upper limits 
to constrain the chemical properties of TN J0924-2201. By comparing the observed emission-line flux ratios with photoionization 
models, we infer that the carbon-to-oxygen relative abundance is already [C/O] > -0.5 at a cosmic age of ~ 1.1 Gyr. This lower limit 
on [C/O] is higher than the ratio expected at the earliest phases of the galaxy chemical evolution, indicating that TN J0924-2201 has 
already experienced significant chemical evolution at z = 5.19. 
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1. Introduction 

The chemical evolution of galaxies is closely related to their past 
star-formation history. Therefore the investigation of the chem- 
ical evolution of galaxies is crucial for understanding the for- 
mation history of galaxies. Observationally, the chemical evo- 
lution of galaxies has been studied by measuring the metallic - 
ity of galaxies at various redshifts. The gas-phase metallicity 
of star-forming galaxies in the local universe can be measured 
fairly easily by means of optical emission-line spectra (e.g., 
[O n]y0727, H/3, [O m] ,15007, Ha, and [N n]^6584; see Nagao et 
al. 2006c and references therein). However, in the early universe 
at z > 1, these diagnostics are more difficult to use because these 
emission lines are often very faint at z > 1 and are redshifted into 
the near-infrared, where sensitive spectroscopic observations are 
more challenging. 

Active galactic nuclei (AGNs), thanks to their high luminos- 
ity, offer an interesting alternative to investigate the metallicity 
in high-z galaxies. The AGN are characterized by a wealth of 
bright emission lines in the wavelength range from rest-frame 
ultraviolet to the infrared, arising from gas clouds photoion- 
ized by their central engines. In particular, the ultraviolet lines, 
e.g., N V/U240, C iV/11549, He il11640, and C m]/U909, are use- 
ful tools to measure the metallicity in high-z AGNs, because 
we can observe these lines in the optical. The metallicity of 
AGNs has been extensively studied (see Hamann & Ferland 
1999 for a review), in particular by focusing on broad-line region 
(BLRs). Nagao et al. (2006a) measured emission-line flux ratios 
of SDSS quasars, e.g., Nv/Civ, (Siiv+Oiv])/Civ, Aim/Civ, 
and Nv/Hen, which are sensitive to metallicity. The authors 
found that there is no redshift evolution of the BLR metallicity 
in the redshift range of 2.0 < z < 4.5. Using near-infrared spec- 
troscopic observations of six luminous quasars at 5.8 < z < 6.3, 



Jiang et al. (2007) found no strong evolution of the BLR metal- 
licity up to z ~ 6. Recently, Juarez et al. (2009) investigated 
the BLR metallicity of a sample of 30 quasars in the redshift 
range of 4.0 < z < 6.4 through the emission-line flux ra- 
tio (Siiv+Oiv])/Civ and found that the metallicity is very high 
even in quasars at z ~ 6. These results indicate that the major 
epoch of chemical evolution in AGNs is at z > 6. However, 
AGN broad lines originate in a very small region in galactic 
nuclei {Rblr < 1 pc; e.g., Suganuma et al. 2006), which may 
have evolved more rapidly and may be not representative of the 
metallicity in their host galaxies. Moreover, most studies have 
reported that the BLR metallicity is significantly higher than the 
solar value (Zblr > Z Q ; e.g., Hamann & Ferland 1992; Dietrich 
et al. 2003; Nagao et al. 2006a; Jiang et al. 2007), reaching as 
much as Zblr ~ 15Z in the most extreme cases (Baldwin et al. 
2003b; Bentz et al. 2004). These extremely high metallicities are 
very hard to reconcile with galaxy chemical evolutionary models 
(e.g., Hamann & Ferland 1993; Ballero et al. 2008) and therefore 
probably reflect the peculiar evolution of the nuclear region. 

To investigate the metallicity on AGN galactic scales, some 
studies have focussed on the narrow line region (NLR). In con- 
trast to the BLRs, the typical size of the NLRs is compara- 
ble to the size of the host galaxies (Rnlr ~ 10 pc; e.g., 
Bennert et al. 2006a, 2006b). The mass of the NLRs is about 
A^nlr ~ 1O 4 ~ 8 M , which is much higher than that of BLRs 
(M BL r ~ 1O 2 ~ 4 M ; see, e.g., Baldwin et al. 2003a). However, 
the metallicity of the NLR in high-z AGN can be constrained 
only in type-2 AGNs, whose BLR and strong continuum emis- 
sion are obscured and therefore allow the detection and careful 
measurement of the narrow lines. Because there are only a few 
optically-selected type-2 quasars discovered at z > 1, many stud- 
ies of the NLR metallicity at z > 1 have focussed on high-z ra- 
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dio galaxies (HzRGs; e.g., De Breuck et al. 2000; Iwamuro et al. 
2003; Nagao et al. 2006b; Humphrey et al. 2008; Matsuoka et al. 
2009). By studying the emission-line flux ratios of N v/C iv and 
Nv/Hen, De Breuck et al. (2000) found that the typical NLR 
metallicity of HzRGs is roughly O.4Z < Znlr < 3.0Z o . On 
the contrary, Nagao et al. (2006b) reported that HzRGs do not 
show any redshift evolution of the NLR metallicity in the red- 
shift range 1.2 < z < 3.8 when exploiting a metallicity diagnos- 
tic diagram involving C iv, He n, and C m], which allows one to 
separate the metallicity and ionization parameter dependencies. 
Matsuoka et al. (2009) confirmed the absence of any significant 
metallicity evolution of the NLRs in HzRGs up to z ~ 4 by using 
a sample of 57 HzRGs at 1 < z < 4. This result suggests that the 
main epoch of major evolution of the NLR metallicity occurred 
at even higher redshifts; it is therefore crucial to investigate the 
metallicity evolution in HzRGs at z > 4. 

In this paper we focus on the most distant radio galaxy, 
TN J0924-2201 at z = 5.19, discovered by Van Breugel et al. 
(1999). TN J0924-2201 is currently the only target that allows 
us to investigate the NLR metallicity at z > 5. In previous stud- 
ies, Venemans et al. (2004) found that TN J0924-2201 is lo- 
cated in a protocluster environment: the density of Lya emitters 
around this object is comparable to the density in protoclusters 
found around HzRGs at z < 4 (see also, Overzier et al. 2006). 
By analyzing the Spitzer data, Seymour et al. (2007) derived the 
stellar mass of TN J0924-2201, i.e., 1O 1IIO M , comparable with 
the stellar masses of HzRGs at 1 < z < 4. Klamer et al. (2005) 
reported that TN J0924-2201 may be a young forming galaxy 
based on its low dust-to-gas ratio (inferred from its Lfir/^co ra- 
tio). This suggests that TN J0924-2201 might be undergoing a 
significant chemical evolution, which can be verified by studying 
the metallicity of its NLR. However, since previous optical spec- 
troscopic observations of this radio galaxy have identified only 
the Lya emission line (i.e., Van Breugel et al. 1999; Venemans et 
al. 2004), we have obtained a very deep optical spectrum of TN 
J0924-2201 with FOCAS at the Subaru Telescope aimed at de- 
tecting additional narrow lines. In the rest of the paper, we adopt 
a concordance cosmology with (Qm, ^a) = (0.27, 0.73) and Ho 
= 71 kms -1 Mpc~'. 



2. Observation and data reduction 

We observed TN J0924-2201 with FOCAS (the Faint Object 
Camera And Spectrograph; Kashikawa et al. 2002) at the Subaru 
Telescope on 2009 Feb. 1 . The observation was performed with 
the 300R dispersion element and the S058 filter to cover the 
5800A < A ohs < 10400A, which cover Lya, N v, C iv, and He n at 
the redshift of the source. We adopted an on-chip binning mode 
of 3 x 1 with a spatial scale along the slit of ~ 07 104 pixeL 1 , be- 
fore pixel binning, and a spectral dispersion of ~ 1 .34A pixeL 1 . 
With a slit width of 078, the resulting spectral resolution was 
8.6A, measured with the observed sky lines. The typical seeing 
was ~ 077. Individual exposure times are 1200s and the total 
integration time is 18000s. 

Standard data reduction procedures were performed by us- 
ing the available IRAF tasks. A bias was subtracted by using 
an averaged bias image, and flat-fielding was performed us- 
ing dome-flat images. Cosmic -ray events were then removed 
by using the lacos_spec task (van Dokkum 2001). The wave- 
length calibration was completed by using sky lines, the sky sub- 
traction by nodding the object along the slit and refined using 
the background task. After sky subtraction, we flux-calibrated 
the spectrum by using observations of the spectrophotometric 
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Fig. 1. Final spectrum of TN J0924-220 1 observed with FOCAS 
after adopting 4 pixel binning in wavelength direction (the mid- 
dle panel). The cyan-shaded regions denote the strong sky-line 
wavelengths. The top panel shows the two-dimensional spec- 
trum. The bottom panel shows the noise spectrum with the same 
spectral binning. 

standard star G191B2B (Oke 1990). After excluding lov/S/N 
frames, we obtained the final stacked two-dimensional spectrum 
(total exposure time is 12000s) and extracted a single spectrum 
of the source using an aperture size of 2718 along the slit. We 
then corrected for the reddening from the Galactic extinction of 
E(B - V) - 0.057 determined from the dust maps of Schlegel et 
al. (1998). 

3. Results 

The final reduced spectrum of TN J0924-2201 is shown in 
Figure 1 . In addition to Lya, the C iv emission is detected with 
high statistical significance (S /N ~ 10). This is the first de- 
tection of the atomic emission line from metals in NLR gas 
clouds at z > 5, which allows us to investigate the chemical 
properties on galactic scales in the early universe. We measured 
emission-line fluxes and profile widths of Lya and C iv by fit- 
ting a Gaussian function using the IRAF task splot (see Table 
1). The line widths after the instrumental broadening correction 
are 1098 and 426 km s , respectively. Note that these widths are 
typical for the NLR emission. For N v and He n lines, 3cr upper 
limits of their fluxes are given by assuming the same line width 
of Civ. 

We calculated the redshifts from the observed wavelengths 
of Lya and Civ. From the Lya line, the derived redshift is 
z = 5.195, which is intermediate between the two previously 
measured redshifts of z - 5.192 (Van Breugel et al. 1999) and 
of z = 5.1989 ± 0.0006 (Venemans et al. 2004). We also derived 
the redshift from Civ line, z - 5.184, which is lower than that 
inferred from the Lya line. The Lya line profile is likely affected 
by absorption from the intergalactic medium (IGM), resulting in 
a artificially higher redshift. Therefore, the redshift derived from 
Civ line (z = 5.184) is more reliable than that from Lya. Note 
that a possible second feature of the Lya line is seen on the red 
wing. This feature may correspond to infalling IGM discussed 
in Dijkstra et al. (2006). 
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Table 1. Emission-line properties of TN J0924-2201. 
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4. Discussion 

4.1. Metallicity 

The detection of the C rv emission allows us to investigate the 
NLR metallicity at z - 5.19. First, the detection of C iv indicates 
that there is plenty of C 3+ ions even at z = 5.19. The CO de- 
tection in Klamer et al. (2005) supports the suggestion that TN 
J0924-2201 had already experienced a significant metal enrich- 
ment in the NLR gas clouds at z > 5.19. 

To investigate a redshift evolution of the metallicity, we com- 
pared some emission-line flux ratios of TN J0924-2201 with 
those of lower-z HzRGs. We compiled emission-line flux data 
of HzRGs from the literature (De Breuck et al. 2000; Matsuoka 
et al. 2009) and plot their Lya/C iv, N v/C iv, and C iv/He n ra- 
tios as a function of redshift (Figure 2). Because the emissivity 
of collisionally-excited emission lines such as C iv increase at 
low metallicity (e.g., Nagao et al. 2006b), the Lya/C iv ratio of 
TN J0924-2201 is expected to be lower than that of lower red- 
shift radio galaxies if there is any significant evolution of the 
NLR metallicity. As shown in Figure 2, the Lya/C iv ratio of TN 
J0924-2201 is somewhat lower than the average value observed 
in radio galaxies at lower redshifts, suggesting some metallic- 
ity evolution. However, the Lya/C iv ratio in TN J0924-2201 is 
still within the scatter observed in lower redshift radio galaxies, 
which prevents us from making any strong claim on the metallic- 
ity evolution based on this line ratio. Moreover, the Lya/C iv flux 
ratio is highly uncertain as a consequence of the IGM absorp- 
tion of Lya, which increases at higher redshift. Another source 
of uncertainty is that the Lya flux may be partly associated with 
star formation in the host galaxy. De Breuck et al. (2000) re- 
ported that Lya is unusually strong in some HzRGs, especially 
at z > 3. Villar-Martin et al. (2007) also confirmed that z > 3 
radio galaxies tend to show higher Lya luminosities than lower- 
Z ones. These results suggest that the radio galaxy phenomenon 
is more often associated with massive starbursts at z > 3 than at 
z < 3. However, TN J0924-2201 does not follow the trend re- 
ported by Villar-Martin et al. (2007). More specifically, the Lya 
luminosity of TN J0924-2201 is 4.8 x 10 43 erg s"\ compara- 
ble with the median of Lya luminosity of HzRGs at z < 3, i.e., 
4.4 x 10 43 erg s" 1 (Villar-Martin et al. 2007). 



The middle panel of Figure 2 shows the emission-line flux 
ratio of N v/C iv as a function of redshift. The N v line is some- 
times regarded as a metallicity indicator for NLRs (e.g., van Ojik 
et al. 1994; De Breuck et al. 2000; Vernet et al. 2001 ; Overzier et 
al. 2001; Humphrey et al. 2008) and therefore it is interesting to 
examine the redshift dependence of the N v/C iv flux ratio as an 
indicator of the chemical evolution of HzRGs. Our upper limit 
on the N v/C iv ratio is below many detections at low redshift, but 
also consistent with several other upper limits. Hence our upper 
limit on N v/C iv cannot really indicate whether or not the metal- 



> 

6 io 



to 



0.1 



0.01 
10 






> 



• This study 

o Matsuoka et al. 2009 

O De Breuck et al. 8000 




O o 



-+- 




Redshift z 

Fig. 2. Flux ratios of Lya/C iv, N v/C iv, and C iv/He n in HzRGs, 
as a function of redshift. TN J0924-2201 is shown with the red- 
filled circle. The blue-open and black-open circles show the pre- 
vious data of HzRGs given in Matsuoka et al. (2009) and De 
Breuck et al. (2000), respectively. The 3<x upper and lower lim- 
its are shown with arrows. 



licity of the NLR of TN J0924-2201 is lower than observed in 
lower redshift galaxies. 

The bottom panel of Figure 2 shows the emission-line flux 
ratio of C iv/He n as a function of redshift. Because the lumi- 
nosity of the C iv line increases at low metallicity, as mentioned 
above, this lower limit indicates that the NLR metallicity of TN 
J0924-2201 at z - 5.19 is comparable to or lower than that at 
z < 5, as in the case of the N v/C iv. 

4.2. Carbon abundance 

In this section we use the observational constraints discussed 
above to investigate the relative carbon abundance in TN 
J0924-2201. The [C/O] abundance ratio is particularly interest- 
ing to investigate galaxy evolution in the early universe, since 
carbon enrichment is delayed compared to a elements because 
a significant fraction of carbon is produced from intermediate- 
mass stars, hence [C/O] can be considered a "clock" of star for- 
mation (e.g., Hamann & Ferland 1993; Matteucci & Padovani 
1993). In particular, a significant decrease in the carbon relative 
abundance should be expected at high redshifts if galaxies are 
indeed in a chemically young phase. 

Motivated by this scenario, we examined the dependence 
of the carbon abundance on the emission-line flux ratios, e.g., 
N v/C iv and C iv/He n, by running photoionization models us- 
ing the Cloudy (version 08.00; Ferland et al. 1998). We used 
the "table AGN" command for the input SED of ionizing pho- 
tons. We then assumed NLR clouds, with hydrogen density in 
the range «h = 10 - 10 4() cm~ 3 , ionization parameters in the 
range U - 10~ 25 - 10~ L5 , and (total) metallicity in the range 
Znlr = 0.5 - 2.0Z o . These parameter ranges are typical of low-z 
AGNs (e.g., Nagao et al. 2001, 2002), though it has not been 
confirmed observationally that these parameter ranges can be 
adopted for NLRs at z > 5. Note that very low or very high 
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Fig. 3. Solid lines show the predicted dependencies of the N v/C iv and C iv/He n flux ratios on the carbon relative abundance [C/O] . 
The panels (a), (b), and (c) show models with (n H , U, Z NhR /Z ) = (10 10 , 1(T 2 - , 1.0), (10 40 , 1(T 2 - , 1.0), and (10 40 , 10" 15 , 0.5), 
respectively. The gray shadows denote the ranges that are consistent withe observational limits on the N v/C iv and C iv/He n flux 
ratios in TN J0924-2201. The allowed [C/O] range is shown as the dotted line and arrow. The minimum carbon relative abundance 
([C/0] m ; n ) for each model is also given in each panel. 



metallicities are not expected for the NLR in TN J0924-2201, 
as discussed in Section 4. 1 . For the chemical composition, we 
assumed that all metals scale by keeping solar ratios except for 
C, He, and N. We then adopted the analytical expressions given 
in Dopita et al. (2000) to obtain the helium and nitrogen relative 
abundances as a function of metallicity. To investigate the depen- 
dence of emission-line flux ratios on the relative carbon abun- 
dance, we examined in our model runs the range -2 < [C/O] 
< 2. Dust grains were not included since we consider relatively 
high-ionization lines (see Nagao et al. 2003). 

Figure 3 shows part of our model results. By comparing our 
observational limits with the model predictions, we infer lower 
limits of the carbon relative abundance ([C/0] m ; n ) for each com- 
bination of parameters used in the Cloudy runs. We obtain that 
carbon in TN J0924-2201 must have been already significantly 
enriched. More specifically, in most cases, with realistic NLR 
physical parameters, the inferred relative carbon abundance is 



weaker star-formation activity and/or heavier IGM absorp- 
tion on Lya. 
3. By comparing observational limits on emission-line flux ra- 
tios (and in particular N v/C iv and C iv/He n) with photoion- 
ization models, we inferred that the carbon is already rel- 
atively high in TN J0924-2201, more specifically [C/O] 
> -0.5, suggesting that TN J0924-2201 is not in an ex- 
tremely young phase, but likely already a few hundred mil- 
lion yrs old. 

Acknowledgements. We thank Yuko Ideue and the Subaru Telescope staff for the 
assisting our FOCAS observation. K.M. acknowledges financial support from the 
Japan Society for the Promotion of Science (JSPS). Data analysis were in part 
carried out on common-use data analysis computer system at the Astronomy 
Data Center, ADC, of the National Astronomical Observatory of Japan. T.N. and 
Y.T. are financially supported by JSPS (Grant Nos. 23244031 and 23654068). 



[C/OL 



-0.5 (e.g., panels a and b in Figure 3). Only in some References 



extreme cases, with low general metal abundance (Z = 0.5Z o ) 
and high ionization parameter (L^ = 10~ L5 ), the inferred lower 
limit may be as low [C/0] m ; n > -0.7. 

The inferred lower limit on the carbon relative abundance 
([C/0] m ; n > -0.5) is higher than some model predictions ex- 
pected for the initial phases of galaxy evolution. For instance, 
Hamann & Ferland (1993) present the evolution of relative ele- 
mental abundance ratios adopting various star-formation histo- 
ries and IMFs, and all of their models predict [C/O] < -0.5 at 
galaxy ages of f age :£ 100 - 300 Myr. Adopting their models, 
our observational limit on [C/O] then suggests that the age of 
TN J0924-2201 must be at least of a few hundred Myr; i.e., ac- 
cording to our analysis TN J0924-2201 is not in the expected 
extremely young phase. 

5. Conclusion 

To investigate the chemical evolution of HzRGs at z > 5, we 
investigated the rest-frame UV spectrum of TN J0924-2201 and 
obtained the following main results. 

1. We detected Civ emission in TN J0924-2201 at z = 5.19; 
this is the first detection of metal emission lines from NLR 
gas clouds at z > 5, suggesting that a significant amount of 
carbon exists even at z > 5 in the early universe. 

2. The measured Lya/Civ flux ratio is lower than that of 
HzRGs at 3 < z < 5; this can possibly be ascribed to 
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